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Abstract 


This  final  report  summarizes  WSU’s  progress  from  4/2/2015  to  09/30/2016  on 
the  project,  “Cross-linked  Metal  Particles  (CLMPs)  for  Low-noise  Bolometer  Ma¬ 
terials”  funded  by  the  US  Army  Research  Office  under  Contract  #  W911NF-15-1- 
0117.  We  have  successfully  carried  out  the  tasks  listed  in  our  proposal  and  gained 
significant  insight  into  electrical  transport  in  these  CLMPs.  We  have  successfully 
synthesized  Au  nanoparticles  with  different  capping  agents,  grown  good  nanoparti¬ 
cle  films,  and  measured  temperature-dependent  resistance  of  the  nanoparticle  films. 
In  addition,  working  with  Prof.  A.  J.  Syllaios  at  the  University  of  North  Texas,  we 
have  started  characterizing  the  noise  performance  of  the  nanoparticle  films.  Our 
results  indicate  that  the  CLMPs  can  simultaneously  have  a  high  temperature  co¬ 
efficient  of  resistivity  and  a  low  noise,  and  therefore  have  a  great  potential  in  the 
microbolometer  technology. 
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1  Introduction 


The  objective  of  this  project  is  to  design,  fabricate,  and  characterize  a  novel  bolometer 
material:  cross-linked  metal  particles  (CLMPs).  In  the  CLMPs,  the  inter-particle  distance 
is  regulated  by  organic  ligands  that  are  chemically  bonded  to  adjacent  metal  particles. 
According  to  our  modeling  results,  the  temperature  coefficient  of  resistance  (TCR)  and 
noise  in  a  CLMP  film  are  determined  by  the  particle  size  and  inter-particle  distance, 
respectively.  Hence  CLMPs  enable  independent  tuning  of  the  TCR  and  1/f  noise  and 
may  revolutionize  the  search  for  appropriate  microbolometer  materials. 

This  project  started  from  the  results  of  an  ARO-funded  project  when  the  PI  was 
employed  at  SRI  International.  Those  previous  results  indicate  that  while  CLMP  films 
are  very  promising  for  low-noise  microbolometers,  growing  thick  CLMP  films  cross-linked 
with  short  ligands  appears  very  challenging.  To  overcome  this  technical  challenge,  at  the 
beginning  of  the  project,  we  designed  several  strategies  and  different  approaches  to  ensure 
a  successful  growth  of  high-quality  CLMP  films.  Using  these  strategies  and  approaches, 
we  have  succeeded  in  synthesizing  Au  NPs  with  several  capping  agents  and  growing  high- 
quality  CLMP  films  cross-linked  with  short  ligands,  as  required  from  modeling  results. 
We  have  also  characterized  the  TCR  of  the  CLMP  films  and,  in  collaboration  with  Prof. 
Syllaios  at  the  University  of  North  Texas,  started  measuring  the  noise  performance.  On 
the  modeling  aspect,  we  have  developed  a  comprehensive  modeling  tool  to  calculate  elec¬ 
trical  transport  and  infrared  absorption  properties  in  CLMP  films  and  to  design  CLMP 
films  for  optimized  performance. 

2  Comprehensive  model  of  electrical  transport  and 
infrared  absorption  of  CLMPs 

The  most  important  material  properties  relevant  to  microbolometer  applications  are 
temperature-dependent  electrical  conductivity,  1/f  noise,  and  infrared  absorption.  In 
this  project  we  have  developed  a  comprehensive  code  that  can  model  these  properties 
effectively.  The  code  includes  the  following  key  ingredients. 

2.1  Coulomb-energy  barrier 

When  the  CLMP  film  is  under  bias,  electrical  transport  occurs  via  electron  tunneling 
between  adjacent  metal  nanoparticles  (NPs).  In  this  tunneling  process,  the  two  originally 
neutral  metal  particles  become  charged:  one  with  a  negative  charge  and  the  other  with  a 
positive  charge.  Thus  the  electron  tunneling  must  overcome  the  Coulomb-energy  barrier 
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between  the  charged  NPs,  the  Coulomb  energy  of  charging  the  two  particles  can  be  ex¬ 
pressed  as  e2/C,  where  C  is  capacitance  between  two  adjacent  particles  and  is  evaluated 
by  considering  an  infinite  series  of  image  charges. 

First  we  examine  the  energy  required  to  oppositely  charge  two  metal  NPs  with  an 
infinite  inter-particle  distance.  Since  an  isolated  metal  sphere  can  be  considered  as  a 
capacitor  with  a  capacitance  of  47t coe^a*,  where  a,  is  the  radius  of  the  ith  sphere  and  e& 
the  effective  dielectric  constant  obtained  above,  the  Coulomb  energy  of  charging  the  two 
separated  NPs  is 


2  Ci  2C2 

e2  /  1  1  \ 

87reoe;>Vai  a2/ 


(2.1) 


The  Coulomb  energy  needed  to  oppositely  charge  adjacent  NPs  would  be  smaller  than 
E®  calculated  above  because  of  the  interaction  between  the  NPs.  Using  the  image-charge 
method,  we  found  that  the  Coulomb  energy  of  charging  the  two  separated  NPs  can  be 
expressed  as 

tp  _  _  e2(fcn  +  k22  +  2 kw) 

2  C  2e/,(knk22  —  kf2) 

where 


OO 

k  11  =  Gqa2  sinh  9  |o2  sinh  nO  +  cq  sinh(n  —  1)0|_1,  (2.3) 

n=  1 

00 

k\2  =  —  ^sinh^y^  Isinhn^C1:.  (2.4) 

c  ^ 

n=  1 
00 

k22  =  «io2  sinh  6  |ai  sinh  nd  +  a2  sinh(n  —  l)9\~l .  (2.5) 

n— 1 


Here  cosh  9  =  (c2  —  a2  —  a2)/(2a!a2)with  a\  and  a2  being  the  radii  of  the  two  spheres,  c 
the  center-to-center  distance  between  the  two  spheres,  c  =  a\  +  a2  + 1  with  t  being  the 
edge-to-edge  distance  (the  ligand  length),  and  is  the  dielectric  constant  of  the  organic 
ligand. 

The  corresponding  capacitance  between  the  two  adjacent  particles  is 


c  _  ebjknkn  ~  fci2) 
kn  +  k22  +  2k\2 


(2.6) 
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Figure  1:  Coulomb  energy  of  oppositely  charged  NPs  as  a  function  of  edge-to-edge  distance 
between  the  particles  for  different  particle  sizes. 

Ec  and  C,  which  depend  on  the  radii  of  and  the  distance  between  the  two  adjacent 
particles,  are  important  properties  that  will  affect  electrical  transport  in  CLMPs.  Figure 
1  plots  the  Coulomb  energy  between  two  particles  as  a  function  of  their  edge-to-edge 
distance. 


2.2  Microscopic  polarizability  of  organic  ligands 

The  Coulomb-energy  barrier  calculated  above  depends  on  the  dielectric  constant  of  the 
host  material  —  in  this  case,  the  organic  ligands  connecting  the  NPs.  However,  these 
ligands  are  molecular  monolayers,  which  renders  the  usual  macroscopic  dielectric  constant 
of  these  organic  materials,  an  ensemble  averaged  property,  inapplicable  to  the  current 
situation.  Instead,  one  must  use  microscopic  polarizability  of  individual  organic  ligands 
to  calculate  the  Coulomb  energy  between  the  charged  NPs. 

We  have  employed  the  Gaussian09  program  suite  to  carry  out  the  polarizability  cal¬ 
culations  for  representative  organic  ligands  that  have  been  used  in  literature  as  cross-link 
agents  between  metal  NPs.  Table  I  lists  the  obtained  polarizability  p  of  these  organic 
molecules  and  their  molecular  volume  v. 

For  molecules  with  a  negligible  permanent  dipole  moment,  the  effective  dielectric  con¬ 
stant  can  be  calculated  from  the  Clausius-Mossotti  equation, 


eb-l  47t 

- =  — p. 

eb  +  2  3v  F 


(2.7) 
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Tabic  1:  Volume  (v),  polarizability  (p),  and  effective  dielectric  constant  (e&)  of  cross-link 
organic  molecules,  ro  is  the  Bohr  radius. 


Molecule 

0  co 

«  (ro) 

Cb 

Ci0N8H6S4O2 

2760.955 

4.846e+02 

9.34 

C12IU0S3 

1942.358 

3.057e+02 

6.08 

C24H18S6 

4285.346 

6.925e+02 

7.29 

C28N2H2O 

3459.573 

8.448e+02 

4.53 

C6H6S2 

1280.082 

1.678e+02 

4.65 

C8IU0S2 

1499.479 

1.922e+02 

4.48 

PPh3 

2433.059 

2.544e+02 

3.34 

2n  -1 

Figure  2:  Polarization  as  a  function  of  length  of  molecules. 


For  molecules  with  a  sizable  permanent  dipole  moment  p,  the  Lorentz-Debye  equation 
should  be  used, 


c±zl  =  il(p+J^) 

eb  +  2  3vV  3  kBT)' 


(2.8) 


where  T  is  the  temperature  and  fc#  the  Boltzmann  constant.  Using  these  equations,  we 
calculate  the  effective  dielectric  constant  e&  for  the  organic  molecules  and  list  them  in 
Table  I.  The  obtained  e&  will  be  used  to  evaluate  the  Coulomb-energy  barrier  of  electron 
tunneling  between  adjacent  NPs. 

We  emphasize  that  the  obtained  e&  would  in  general  be  different  from  the  dielectric 
constant  of  a  thick  him  of  these  molecules.  In  fact,  the  polarizability  is  not  proportional 
to  the  system  size,  as  shown  in  Figure  2,  which  plots  the  polarizabilities  of  polyyne 
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and  polyene  molecules,  both  chain-like,  versus  the  molecular  length.  It  is  clear  that  the 
polarizability  depends  on  the  molecular  length  nonlinearly. 

2.3  Electron  tunneling  between  two  adjacent  particles 

To  study  electrical  transport,  we  consider  the  basic  continuity  equation  (charge  conserva¬ 
tion)  at  a  given  metal  particle, 


- u  <2-9> 

3 

where  nt  is  the  electron  density  in  the  ith  particle,  and  ItJ  is  the  electrical  current  flowing 
from  particle  %  to  particle  j. 

When  a  CLMP  film  is  under  a  bias,  individual  particles  have  a  different  potential,  Vt 
for  particle  i.  Then  the  electron  tunneling  current  from  occupied  states  in  particle  i  to 
unoccupied  states  in  particle  j  is 

hj  =  f  dEpe~2xtii  f(E)[l  -  f(E  +  q(Vj  -  K)  ~  E^)\,  (2.10) 

where  f(E)  =  [1  +  eE/fcsT]-1  is  the  Fcrmi-Dirac  distribution,  p  is  normalization  constant, 
tij  is  the  edge-to-edge  separation  between  particles  i  and  j,  E^  is  the  Coulomb  energy  of 
charging  particles  i  and  j  obtained  earlier,  and  x  —  [2 m(<fr—E)]1//2/ h  is  the  tunneling  factor 
obtained  from  the  WKB  approximation  with  (j)  and  m  being  the  metal  work  function  and 
the  electron  mass.  Here  we  set  the  electrochemical  potential  (Fermi  energy)  of  the  particle 
i  to  be  zero.  This  expression  describes  two  competing  factors  that  influence  the  electron 
tunneling:  One  is  the  potential  difference  between  the  two  particles  in  the  presence  a  bias, 
which  facilitates  electron  tunneling;  the  other  is  the  Coulomb  energy  of  charging  the  two 
particles  (because  after  the  tunneling,  particle  j  becomes  negatively  charged  and  particle 
i  becomes  positively  charged),  which  hinders  electron  tunneling. 

Conversely,  the  electron  tunneling  current  from  particle  j  to  particle  i  is 

hi  =  f  dEpe~2xt':i  f  {E  +  q(Vj  -  E))[l  -  f(E  -  E^)}.  (2.11) 

The  basic  tunneling  probability  is  evaluated  at  finite  temperature  and  at  arbitrary 
voltage  difference. 
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2.4  Electrical  transport  in  a  network  of  particles 


To  calculate  the  I-V  characteristics  of  the  CLMP  him,  we  need  to  consider  a  3d  network 
of  metal  particles  and  solve  the  steady-state  equations, 


where 


1 

;o> 

UP 

O 

(2.12) 

13  p  i  _  emvj-v^)-E9. \  > 

(2.13) 

l  _  -2 xtii  q(V3  +  Eij 

31  p  e/3  kw-w+effi-i' 

(2.14) 

These  equations  are  highly  nonlinear  and  the  number  of  equations  is  enormous  (we  are 
currently  doing  32x32x32). 

Under  a  low  bias,  the  expressions  of  IVJ  and  Iji  can  be  linearized.  The  obtained  linear 
equations  can  be  efficiently  solved  by  using  the  sparse-matrix  technique. 

Obtaining  a  stable  solutions  for  these  nonlinear  equations  is  very  challenging.  We  have 
developed  a  numerical  algorithm,  in  which  we  gradually  increase  the  applied  voltage 


vf  =  v°  +  sv,. 


(2.15) 


where  V)°  is  the  solution  to  the  equations  for  an  applied  bias  of  V°,  and  V)  is  the  solution 
to  be  found  when  the  bias  increases  by  SV,  V  =  V°  +  SV.  We  make  qSV  <C  /cbT,  and 
thus  SVi  is  small  as  compared  with  Id0 .  We  linearize  the  equations  around  ld°,  and  the 
above  nonlinear  equations  become  linear  equations  of  SVi 


£ 

3 


pe~2xUj  (SVj  -  SV) 


e-P[q(y^-vP)-E9.] 


e/%(U0-U°)+Ei,] 


q2(3[q(V f  -  E°)  -  q2(3[q( V?  -  Id0)  +  Eci3]ep[q{y? ~V^+El 


[l-e 


J2pe 


-2  XU 


-mv?-vr>)-E^  2 
q(V?  -  V?)  -  E\ 


-mvq-vp)-E9.] 


L1  -  e 


+  pe 


[eP\q^-v?)+E9. )  _  1p 

1(K  -  v«)  +  Et, 


MVf-vn+E: 


J3_ 
ij\  —  1 


(2.16) 


which  can  be  solved  by  using  the  sparse  matrix  techniques  and  the  conjugated  gradient 
method. 

Figures  3  and  4  plot  current  as  a  function  of  the  inverse  of  applied  bias  for  CLMP 
films  at  different  temperatures.  It  is  seen  that  the  I-V  curves  are  highly  nonlinear  as 
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Figure  3:  Current  vs  inverse  of  bias  in  CLMPs  for  e&  =  2.  Black,  red,  and  green  curves 
correspond  to  T  =  300,  200,  and  100  K,  respectively.  The  parameters  are  a  =  2  nm  and 
t  =  0.5  nm. 


the  applied  bias  increases.  There  is  a  clear  crossover,  above  which  the  I-V  characteristics 
become  independent  of  temperature.  This  crossover  occurs  when  the  electron  energy 
gain  due  to  the  voltage  drop  overcomes  the  Coulomb  energy  of  charging  two  adjacent 
particles.  Thus  the  Coulomb  energy  can  be  obtained  from  the  voltage  dependence  at  a 
fixed  temperature  and  from  the  temperature  dependence  under  a  small  voltage.  In  the 
high-bias  regime,  the  logarithm  of  current  depends  linearly  on  1/C,  suggesting  that 

/~exp(— C/C),  (2.17) 

where  the  characteristic  bias  Vc  is  the  bias  at  which  the  Coulomb  energy  (£/))  is  the  same 
as  energy  gain  due  to  the  potential  difference,  +  cij  +  t)),  with  L  being  the  length 

between  the  two  electrodes. 

Figures  3  and  4  also  show  that  strong  temperature  dependence  of  electrical  transport  is 
at  low  bias,  and  that  to  achieve  high  temperature  sensitivity,  a  large  (£)/■)  is  needed.  Thus 
the  Coulomb  energy  can  be  obtained  from  the  voltage  dependence  at  a  fixed  temperature 
and  from  the  temperature  dependence  under  a  small  voltage. 

2.5  AC  transport 

The  AC  transport  introduces  another  time  scale,  i.e. ,  the  time  period  the  AC  voltage. 
In  addition,  at  a  finite  frequency,  both  electrical  tunneling  and  displacement  current  in  a 
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Figure  4:  Current  vs  inverse  of  bias  in  CLMPs  for  eb  =  6.  Black,  red,  and  green  curves 
correspond  to  T  =  300,  200,  and  100  K,  respectively.  The  parameters  are  a  =  2  nm  and 
t  =  0.5  nm. 


capacitor  between  any  two  adjacent  particles  contribute  to  electrical  transport.  Thus  the 
AC  transport  would  also  reveal  information  on  the  capacitance  between  adjacent  particles, 
which  depends  on  the  particle  size  and  the  inter-particle  distance. 

To  study  the  AC  transport,  one  need  to  solve  time-dependent  equations. 

9^  =  £«<  -  hi)-  (2.18) 

3 

At  low  bias,  |  q(Vj  —  V.)  |  -C  ksT,  the  AC  transport  problem 
network  problem.  For  an  AC  bias  with  frequency  of  oj,  V (t)  = 
reduce  to 

0  =  £(U  -  Vi)Zy‘, 

3 

with  an  impedance  between  each  pair  of  adjacent  particles, 

=  RiJ  +  iwCij, 

c  _  ebkuk22  ~  kf2 
^  ku  +  k22  +  2k12 
2 

—  rip  —  p~ XUj _ 

O  P  i  _  e/3£% 


(2.20) 

(2.21) 

(2.22) 


is  reduced  to  an  impedance 
Ve~lut,  the  above  equations 

(2.19) 
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In  deriving  Eq.  (2.20),  we  have  neglected  the  Fermi  energy  change  in  each  metal  particle 
of  the  CLMP  him  under  a  small  AC  voltage. 

To  solve  these  complex-coefficient  linear  equations,  we  first  treat  the  real  and  imaginary 
parts  of  complex  V)  as  two  independent  real  variables  and  rewrite  these  equations  as  real- 
coefficient  linear  equations.  The  matrix  size  changes  from  NxN  to  2N x 2N.  We  fold  that 
the  conjugated  gradient  approach  is  not  numerically  reliable.  So  we  have  employed  the 
generalized  minimal  residual  algorithm,  together  with  sparse-matrix  technique  to  solve 
these  equations  and  ford  complex  V).  Then  we  compute  the  total  current  /  and  obtain 
the  AC  conductance  cr(co)  via  cr(to)  =  I  (to) /V  (to). 

Figure  5  plots  both  the  real  part  and  the  modulus  of  cr  as  a  function  of  frequency  c o. 
We  see  that  |cr(u;)|  depends  only  weakly  on  the  frequency  at  low  frequencies,  then  takes 
off  at  a  critical  frequency,  increasing  linearly  with  the  frequency.  The  underlying  physics 
is  that  at  low  frequencies,  the  displacement  current  through  the  capacitance  is  negligible 
as  compared  to  the  tunneling  current,  electrical  transport  is  due  mainly  to  tunneling, 
which  is  independent  of  frequency.  As  the  frequency  increases,  the  displacement  current 
becomes  more  important,  and  at  the  critical  point,  is  comparable  to  the  tunneling  current. 
After  that,  the  displacement  current  is  dominant.  In  the  figure,  the  slope  of  |(j(a;)|  at 
large  to  is  determined  by  the  effective  capacitance  of  the  system.  Thus  from  impedance 
measurement  of  an  CLMP  him,  we  should  be  able  to  extract  the  information  of  the  average 
capacitance  and  tunneling  current  between  adjacent  particles.  Since  they  both  depend 
on  the  particle  size  and  the  inter-particle  distance,  the  AC  transport  measurement  could 
serve  as  a  diagnostic  tool  to  independently  determine  the  particle  size  and  volume  loading 
density  in  a  CLMP  him. 


2.6  Noise 


Noise  is  critical  to  the  performance  of  a  microbolometer.  In  fact,  the  sensitivity  of  existing 
microbolometers  based  on  VCL  and  a-Si  is  limited  by  the  1/f  noise  in  these  systems. 
In  general,  the  noise  originates  from  the  diverse  timescales  associated  with  the  electron 
tunneling  between  adjacent  NPs. 

The  model  also  systematically  calculates  the  noise  level  by  solving  the  time-dependent 
master  equation,  in  the  present  of  a  temporal  fluctuation,  r]i(t), 


where 


dfi/ii 

dt 


Y^lwij(SVi  -  <W1 

3 


W. 


1-3 


7 


„2p^2 

e  1  a  c-2Ttc-E1,?  /kBT 

Ant 


(2.23) 


(2.24) 
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Figure  5:  AC  conductance  of  CLMP  film  as  a  function  of  frequency.  Circles  and  squares 
correspond  to  Recr(o;)  and  |er(u;)|,  respectively. 


The  master  equation  can  be  written  as  a  matrix  form, 

dS^i 


dt 


^  1  AjjSfij  rji , 


(2.25) 


where  the  matrix  A  is 


Aij  ^  ^  Wik.dijj  T 


Wij. 


(2.26) 


Denoting  the  kth  eigenvalue  of  A  as  A  &  and  its  associated  eigenvector  as  the 
spectral  density  of  fluctuations,  according  to  the  Wiener-Khintchinc  theorem,  is 


Pk{u)  =  2  £k{u)£k(-u)  = 


2  Dk 


CO2  +  XI 


(2.27) 


where  D *.  is  the  averaged  intensity  of  the  random  pulses,  and  can  be  determined  by  the 
equipartition  law,  that  the  chemical-potential  (voltage)  fluctuation  in  a  circuit  is 


\a$  =  \kBT= fpt(u)<L,=  !±, 


(2.28) 


where  C  is  the  average  capacity  of  the  circuit.  The  total  spectral  density,  including  all 
the  modes,  is  then 

(2.29) 
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Figure  6:  Distribution  of  logarithm  of  eigenvalue  (left  panel)  and  noise  spectral  density 
for  CLMPs  with  different  size  fluctuation  density. 


Thus  to  calculate  the  noise  spectrum,  we  need  to  find  all  the  eigenvalues  of  the  matrix 
A.  It  is  easy  to  see  that  one  eigenvalue  is  A0  =  0,  which  corresponds  to  the  equilibrium 
situation.  All  other  eigenvalues  must  be  negative,  A*,  <  0,  because  any  fluctuation  will 
decay  over  time.  We  have  considered  a  large  CLMP  system  consisting  of  32  x  16  x 
16  metal  NPs,  and  computed  all  eigenvalues  A&  as  well  the  noise  spectral  density  for 
different  particle  size  fluctuations,  as  shown  in  Figure  6.  We  see  that  1/Afc,  which  describes 
microscopic  time  scales,  has  a  very  narrow  distribution  as  compared  to  disordered  VOx. 
Hence  the  noise  level  in  CLMPs  is  much  lower  than  that  in  disordered  materials  like 
VOx.  More  important,  the  noise  level  becomes  saturated  at  low  frequencies  and  does  not 
increase  indefinitely  as  the  frequency  is  lowered.  Thus  the  CLMPs  would  have  a  very  low 
noise  level. 

2.7  Infrared  absorption 

Since  a  microbolometer  is  used  to  detect  infrared  radiation,  it  is  important  to  understand 
the  IR  response  of  the  CLMP  films  between  8  to  14  microns.  When  the  active  bolometer 
material  is  a  good  IR  absorber,  the  microbolometer  architecture  can  be  greatly  simplified. 
Our  model  also  computes  the  IR  response  of  CLMP  films,  which  can  be  considered  as 
metal  NPs  embedded  in  an  organic  host.  The  effective  dielectric  response  of  this  composite 
system  usually  can  be  computed  from  the  well-known  Maxwell- Garnett  theory  and  the 
Mie  solution. 
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Table  2:  n  and  k  in  the  infrared  region  of  An  bulk 


wavelength  (pm) 

n 

k 

7.999 

8.5 

46.4 

8.266 

9.016 

47.6 

8.551 

9.582 

48.8 

8.856 

10.21 

50.2 

9.184 

10.84 

51.6 

9.537 

11.51 

53.1 

9.919 

12.24 

54.7 

eeff(tA)  -  €i(u)  (l  +  ),  (2.30) 

£p(u})-£  l(w)  •> 

where  e±  is  the  dielectric  response  of  the  host  (organic)  material,  /  is  the  volume  fraction 
of  NPs  in  the  composite, 


f  47 t  a3 

*  ~  T  (2 a  +  f)3’ 

and  ep  is  related  to  the  dielectric  function  of  metal  62  via 


(2.31) 


^  =  F(x) 
£2 


2  (sin  a;  —  xcosx) 
x2  —  1)  sin  x  +  x  cos  x  ’ 


(2.32) 


with  x  =  ka^/e^JF  and  k  being  the  wave  vector. 

We  use  the  measured  dielectric  function  of  bulk  Au  to  describe  e2.  In  literature,  the 
available  n  and  k  in  the  infrared  region  for  Au  are  listed  in  Table  2.  We  interpolate  these 
measured  values  to  obtain  n  and  k  for  other  wavelengths. 

Special  caution  must  be  taken  when  describing  the  dielectric  function  of  the  host 
(organic)  material.  Usually,  an  insulating  bulk  organic  material  does  not  strongly  absorb 
the  IR  radiation.  In  CLMPs,  however,  the  organic  ligand  is  not  an  electric  insulator,  for 
electrons  can  tunnel  between  adjacent  metal  NPs.  Thus  the  organic  material  in  CLMPs 
do  have  an  effective  conductivity.  In  addition,  as  discussed  earlier,  there  is  a  capacitance 
C  between  each  pair  of  adjacent  NPs,  which  also  facilitates  electrical  conduction  at  finite 
frequencies.  Hence  the  effective  resistivity  of  the  organic  ligand  between  the  two  adjacent 
NPs  at  frequency  u  is 

2  at  (  ^  t  \  -1  , 

pe  =  ^Tv1  +  lujCPt^)  >  (2-33) 
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Figure  7:  IR  absorption  spectrum  of  CLMPs.  The  parameters  are  a  =  2  nm  and  t  =  0.5 
nm. 

where  A  is  the  effective  tunneling  area  with  the  non-planar  nature  of  spheres  included, 
A  =  Tra/r ,  and  pt  =  1/cr  with  a  from  Eq.  (2.22). 

The  equivalent  dielectric  function  of  the  organic  ligands  in  CLMP  films  is  then 


(2.34) 


Substituting  Eq.  (2.34)  into  Eq.  (2.30),  we  obtain  the  effective  dielectric  function  for 
the  composite,  the  CLMP  films.  The  absorption  spectrum  can  be  computed  from 


(2.35) 


The  calculated  IR  absorption  spectrum  is  illustrated  in  Figure  7.  We  see  that  the  IR 
absorption  is  greatly  enhanced  when  the  tunneling  conduction  is  included.  Without  the 
tunneling  contribution,  as  illustrated  by  the  green  line,  the  IR  absorption  would  be  very 
weak.  Our  results  suggest  a  strong  IR  absorption  in  CLMP  films  where  electron  tunneling 
between  adjacent  NPs  is  significant. 

3  Synthesis  of  NPs  and  growth  of  NP  films 

In  the  proposed  CLMP  structure,  metal  particles  are  separated  by  organic  ligands  that 
form  chemical  bonds  to  the  metal  particles  at  the  both  ends  (cross-linking).  We  have  tried 
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Figure  8:  Chemical  structures  of:  (a)  ethanethiol,  (b)  propanethiol,  (c)  butanethiol,  (d) 
pentanethiol,  and  (e)  hexanethiol. 


s 
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Figure  9:  UV/Vis  absorption  spectrum  of  Au-Cg  nanoparticles  in  hexanes.  Inset  shows 
the  Au  NP  dispersion,  (b)  The  corresponding  FTIR  spectrum  in  a  KBr  pellet. 


many  different  synthesis  and  film-growth  routes.  Here  we  summarize  the  approaches  that 
worked.  For  an  easy  reference,  we  display  in  Figure  8  the  molecular  structures  of  the 
capping  agents  used  in  this  report. 

3.1  Synthesis  of  Au  NPs  capped  with  1-hexanethiol 

We  synthesized  Au  NPs  from  HAUCI4.3H2O  and  capped  these  NPs  with  1-hexanethiol. 
The  obtained  capped  NPs,  named  as  Au-C6,  are  purified  by  multiple  precipitation/re¬ 
dispersion  cycles  with  toluene  as  the  solvent  and  ethanol  as  the  precipitating  medium. 

Figure  9a  shows  the  UV/Vis  spectrum  of  the  1-hexanethiol  capped  Au  NPs  in  hexanes. 
A  weak  hump,  caused  by  the  plasmon  resonance  of  the  Au  NPs,  is  seen  between  400  and 
600  nm.  The  weak  plasmon  resonance  is  consistent  with  the  small  diameter  of  NPs.  Figure 
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9b  shows  the  corresponding  FTIR  spectrum.  The  observation  of  the  prominent  peak  at 
1462  cm-1  confirms  the  presence  of  1-hexanethiol.  The  absence  of  the  S-H  vibrational 
peak  between  2600  and  2550  cm-1  indicates  that  1-hexanethiols  are  adsorbed  onto  the 
Au  NPs  surface. 

3.2  Synthesis  of  Au  NPs  capped  with  ethanethiol 

From  our  modeling  results,  in  the  Au  NP  film  cross-linked  with  1,9-nonanedithiol,  the 
distance  between  Au  particles,  which  is  the  molecular  length  of  1,9-nonanedithiol,  is  too 
long  for  the  film  to  be  a  semiconductor.  For  the  bolometer  applications,  we  need  to  cross 
link  the  NPs  with  much  shorter  ligands,  like  benzene-1, 4-dithiol,  which  requires  Au  NPs 
capped  with  a  short  agent,  such  as  ethanethiol,  be  synthesized. 

We  tried  to  synthesize  the  ethanethiol  functionalized  Au  NPs,  named  as  AU-C2,  in  the 
same  way  as  Au-C6  described  above,  except  that  chloroform  is  used  as  the  solvent.  The 
use  of  a  shorter  ligand,  ethanethiol  in  place  of  1-hexanethiol  in  Au-C6,  however,  yielded 
much  larger  sized  NPs,  as  indicated  by  the  more  pronounced  plasmon  absorption  between 
400  and  600  nm. 

Facing  the  upset,  we  tried  another  approach  among  our  developed  strategies  and 
approaches.  That  is  to  first  synthesize  triphenyl  phosphine  (TPP)  capped  NPs,  named 
as  Au-TPP.  The  Au-TPP  then  could  undergo  a  ligand  exchange  with  ethanethiol,  giving 
rise  to  Au-C2. 

The  protocol  for  preparing  Au-TPP  is  done  according  to  the  procedure  reported  by 
Woehrle  [2],  Au-TPP  is  formed  by  addition  of  NaBH4  to  a  solution  of  AuCl(TPP)  in 
ethanol.  The  product  is  isolated  by  precipitation  in  hexanes,  followed  by  filtration.  It  is 
then  rinsed  with  hexanes,  CF^CR/hexanes  (1:1),  and  CH2Cl2/hexanes  (3:1). 

Figure  10  shows  UV/Vis  spectrum  of  the  obtained  Au-TPP  in  a  dilute  CHC13  solution. 
The  weak  hump  between  500  and  550  nm  suggests  a  small  diameter  of  the  Au  NPs.  The 
overall  appearance  of  the  spectrum  is  in  agreement  with  that  reported  in  the  literature 
[2]-  . 

To  make  Au-C2,  the  obtained  Au-TPP  NPs  are  dissolved  in  chloroform  and  then 
ethanethiol  is  added  into  the  solution  at  equal  mass  with  Au-TPP.  The  product  is  isolated 
via  precipitation  in  hexanes.  Figure  10b  compares  the  UV/Vis  spectrum  of  Au-TPP  with 
Au-C'2 •  The  weakening  of  the  peak  at  308  nm  indicates  the  removal  of  the  TPP  group  in 

Au-C2. 
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Figure  10:  (a)  UV/Vis  absorption  spectrum  of  Au-TPP  in  a  dilute  CHC13  solution,  (b) 
Comparison  of  UV/Vis  spectra  of  Au-TPP  and  AU-C2. 

3.3  Growth  of  CLMP  films 

To  ensure  good  film  formation,  the  glass  substrate  was  treated  with  piranha  solution  (1/3 
sulfuric  acid/hydrogen  peroxide  mixture)  to  enrich  the  surface  with  hydroxyl  groups.  It 
is  then  rinsed  with  abundant  deionized  water  and  acetone.  Next,  the  glass  substrate  was 
functionalized  with  n-octyltriethoxysilane. 

Following  the  protocol  described  in  [1],  25  /iL  of  Au-Cg  solution  in  heptane  is  drop 
casted  onto  the  n-octyltriethoxysilane  functionalized  glass  substrate.  The  spin-coated  film 
is  then  cross-linked  by  soaking  in  a  solution  of  1,9-nonanedithiol.  Figure  11  shows  the 
obtained  Au-Cg  films  and  the  small  particle  size  is  corroborated  by  the  UV/Vis  absorption 
spectra. 

To  grow  good  A11-C2  films,  again,  the  treatment  to  the  glass  substrate  is  essential.  We 
then  spin-coated  the  precipitated  Au-TPP  on  the  n-octanetriethoxysilane  functionalized 
glass  substrate.  The  film  was  subsequently  treated  with  benzene  1,4-dithiol,  the  desired 
short  ligand,  for  cross-linking. 

Figure  12  shows  the  UV/Vis  absorption  of  the  cross-linked  film.  The  absorption  peak 
within  500  and  550  nm  is  not  seen,  suggesting  that  the  Au  NPs  diameter  is  small.  The 
film  image,  displayed  in  the  inset,  demonstrates  the  quality  and  uniformity  of  the  film. 
Figure  12b  is  the  corresponding  reflectance  spectrum,  which  can  be  fitted  to  allow  film 
thickness  estimation.  The  estimated  film  thickness  is  around  several  hundred  nanometers. 

To  estimate  the  thickness  of  the  films,  we  used  the  Atomic  Force  Microscopy  (AFM) 
and  measured  the  him  thickness.  The  thickness  range  is  between  100-500  nm. 
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Figure  11:  Pictures  of  Au-C@  film  spin-coated  on  an  n-octyltriethoxysilane  functionalized 
glass  substrate,  subsequently  cross-linked  with  1,  9-nonanedithiol  in  (a)  1-butanol  and 
(b,c)  isopropanol,  (d)  The  corresponding  UV/Vis  absorbance  spectra. 
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Figure  12:  (a)  UV/Vis  spectrum  of  A11-C2  spin  coated  from  a  50mg/mL  solution  in 
CHC13,  subsequently  treated  with  benzene  1,4-dithiol  in  hexane,  (b)  The  corresponding 
reflectance  spectrum. 
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Figure  13:  TEM  images  of  A11-C2  taken  at  (a)  80K  magnification  and  (b)  size  distribution. 

3.4  TEM  Images  and  size  of  Au  NPs 

Having  made  the  films  of  Au-Cg  and  AU-C2,  it  is  critical  to  check  on  the  NPs  size. 
The  three  types  of  Au-NPs  that  had  been  prepared  i.e.,  (1)  Au-C2,  (2)  Au-C6,  and  (3) 
Au-TPP  were  sent  out  for  TEM  imaging  at  Franceschi  Microscopy  &  Imaging  Center 
(WSU-Pullman) . 

3.4.1  Au-C2 

Direct  particle  size  measurements  from  its  TEM  images  shows  that  the  diameter  of  A11-C2 
is  about  5.7  ±  1.2  nm  (Figure  13).  The  AU-C2  NPs  are  obtained  by  performing  ligand 
exchange  of  Au-TPP  with  ethanethiol.  Au-TPP  is  known  to  be  unstable  due  to  the  weak 
interaction  with  the  capping  agent  TPP.  It  is  possible  that  the  NP  size  grows  from  the 
initial  <  5  nm  size  to  about  5.7  ±  1.2  nm  either  during  storage  or  during  the  ligand 
exchange  reaction. 

3.4.2  Au-C6 

Due  to  the  poor  contrast,  the  particle-size  measurement  is  less  straightforward  for  Au-Ce, 
but  from  Figure  14  below  it  is  seen  that  the  diameters  are  less  than  5  nm. 

3.4.3  Au-TPP 

Based  on  limited  number  of  particles,  the  average  particle  size  derived  from  Figure  15 
below  is  for  Au-TPP  is  about  2.4  ±  1.0  nm.  This  size  is  ideal  for  bolometer  applications 
from  our  modeling  results.  The  increase  in  Au  NP  size  from  Au-TPP  to  A11-C2  may 
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Figure  14:  TEM  images  of  Au-C6  obtained  at  240K  magnification. 


take  place  during  either  the  storage  or  the  ligand  exchange  process,  as  mentioned  above. 
Currently  we  are  working  on  the  conditions  that  prevent  such  size  change  from  happening. 


4  Electrical-transport  measurements 

We  have  put  together  a  PID  (Proportional,  Integral,  Derivative)  controlled  heater  strip 
to  enable  temperature  controlled  experiments  between  ambient  temperature  to  about  500 
K. 

We  have  measured  the  I-V  curves  using  pre-patterned  electrodes  with  10  /jrn  inter- 
digitated  spacing  (NanoSPR,  LLC.).  A  pre-patterned  is  dipped  into  a  100  mg/mL  gold 
nanoparticle  solution,  and  the  dried  him  is  measured  using  a  Keithley  2400  source  meter 
(Tektronix,  Inc.)  in  the  absence  of  light  in  a  black  box.  The  temperature  of  the  him 
is  varied  from  room  temperature  and  383  K  using  a  custom-built  heater  system,  which 
consists  of  a  strip  heater,  a  temperature  controller,  and  a  thermocouple.  The  I-V  char¬ 
acteristics  and  the  temperature-dependent  resistance  for  Au-C6,  Au-C4,  and  Au-C2  films 
are  plotted  in  Figures  17-19. 

The  resistance  of  the  Au-C6  him  R(T )  changes  from  3.0  MO  at  293  K  to  2.25  MO  at 
303  K,  which  translates  to  a  TCR  of 

1  d  /? 

TCR  =  =  (3.0  -  2.25)/(3.0  *  10)  =  2.5%.  (4.1) 

R  dT 
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Figure  15:  TEM  images  of  Au-TPP  obtained  at  200K  magnification. 


Figure  16:  Schematic  of  apparatus  used  for  I-V  characteristic  measurements. 
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Figure  17:  I-V  characteristics  for  different  temperatures  (a)  and  device  resistance  as  a 
function  of  temperature  (b)  of  Au-Cg  (rgAR044-23)  film. 
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Figure  18:  I-V  characteristics  for  different  temperatures  (a)  and  device  resistance  as  a 
function  of  temperature  (b)  of  A11-C4  (rgAR069-ll)  film. 
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Figure  19:  I-V  characteristics  for  different  temperatures  (a)  and  device  resistance  as  a 
function  of  temperature  (b)  of  Au-C2  (rgAR047-37a)  film. 
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Figure  20:  Schematic  of  apparatus  used  for  noise  measurements. 


This  TCR  is  comparable  to  the  state-of-the-art  value  of  VO,,..  Other  films  also  exhibit  a 
large  TCR. 


5  Noise  measurements 


A  schematic  of  the  noise  measurements  is  shown  in  Figure  20.  The  measured  noise  is 
defined  as 

Vl  =  V(Rf)2  +  (G14)2  (5.1) 

where  V(Rp)  is  the  (Johnson)  noise  of  the  amplifier  feedback  resistor,  Vn  is  the  sample 
noise  and  G  is  the  amplifier  gain.  The  gain  is  the  ratio  of  the  amplifier  output  signal 
voltage,  Vout,  to  the  input  signal  voltage,  V)n,  which  is  equal  to  the  ratio  of  the  feedback 
resistance  to  the  sample  resistance,  R. 

G  =  Vout/Vin  —  Rf/R-  (5-2) 


The  system  noise  is  also  characterized  by  using  metal  him  resistors  whose  resistances 
are  comparable  to  the  devices  tested.  Using  the  metal  him  resistor  noise  results,  the  low 
frequency  noise  of  the  cross-linked  metal  particles  him  noise  and  system  noise  can  be 
quantihed. 

In  general,  the  total  noise  measured  is  comprised  of  1//  noise,  white  noise,  and 
Generation-Recombination  (G-R)  noise,  which  can  be  expressed  quantitatively  as 


Vnoise  _  //  ,  Vwhite  ,  ^G-R 

A/  A/  +  A/  +  A/  • 

The  1//  noise  can  be  described  using  the  Hooge  model,  with 

v_m  a»  Yki  -  a 

A/  pV  p  p ' 


(5.3) 


(5.4) 
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where  p  is  the  carrier  density,  V  is  the  device  volume,  Vbu1s  is  the  device  bias,  and  7  is 
the  frequency  exponent,  typically  near  a  value  of  1  (in  the  range  from  0.9  to  1.1).  The 
normalized  Hooge  parameter  is  olu/p  with  an  being  the  Hooge  coefficient.  For  the  Hooge 
model,  \/~B  oc  Vbms  and  B /V^ias  oc  1/V  for  fixed  bias  voltages. 

The  G-R  noise  for  a  single  trap  is  described  by 

Yg-r  _  ATVgias 
A/  1  +  (2t tt/)2’ 


where  r  is  the  trap  time  constant  and  A  =  nt/Vp2  with  nt  being  the  trap  density.  The 
noise  waveform  for  trapping  can  be  analyzed  by 


r  i^noise  jj white  tt  ,  f  At 

VLs  ~PV  1  +  (2tt/)2’ 


which  can  be  normalized  to  the  device  volume  with 


WN  = 


f 


(K 


2 .  —  v2 

noise  y  white 


v2 

v  Bias 


v 


«  f(nt/p2)T 
p  l  +  (27 r/)2‘ 


(5.6) 


(5.7) 


Here  it  is  assumed  that  the  exponent  of  the  1  //  noise  frequency  is  equal  to  1.  For  the 
case  of  a  single  trap,  a  plot  of  of  Wjy  (determined  from  the  measured  noise  data)  yields 
a  horizontal  line  (of  value  aji/p  )  superimposed  with  a  Lorentzian  (trap).  Figure  23 
shows  an  example  of  Wn  plot  as  a  function  of  frequency  for  noise  data  taken  at  different 
temperatures.  For  temperatures  of  360  K  and  380  K,  the  curves  have  no  slope,  indicative 
that  the  measured  low  frequency  noise  is  1  //  noise  only.  The  values  of  these  curves  are 
equal  to  an/p •  For  temperatures  400  K  and  above,  the  curves  at  low  frequency  (less  than 
10  Hz)  increase  with  increasing  temperature. 

Experimental  room  temperature  noise  measurements  are  made  in  a  shielded  and  light 
tight  enclosure.  For  temperature  dependent  noise  measurements,  the  him  is  mounted 
within  a  cryostat.  All  measurements  are  made  with  the  films  in  dark  conditions.  For 
noise  measurements,  a  Keithley  428  source  is  used  to  bias  the  device  and  as  a  low  noise 
current  amplifier.  The  gain  is  adjusted  by  use  of  a  variable  feedback  resistor  within  the 
preamplifier.  The  amplified  signal  is  converted  to  a  noise  voltage  -  frequency  spectrum 
by  a  HP  35670A  dynamic  signal  analyzer  (DSA).  The  data  is  then  stored  on  computer 
for  analysis. 

Noise  measurements  made  on  two  different  groups  of  cross-linked  metal  particle  films 
Au-C6  and  Au-C6  (TPP2)  deposited  on  interdigitated  electrodes  as  shown  in  Figure  22. 
The  noise  measurements  were  made  at  Room  temperature.  The  noise  data  for  the  films 
with  metal  particles  rgAR081-la  Au-C6  is  shown  in  Figure  23.  For  frequencies  greater 
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Figure  21:  Normalized  waveforms  analyzed  for  trapping  at  various  temperatures. 


Figure  22:  Sample  film  on  a  NanoSPR  interdigitated  electrode. 
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Figure  23:  Noise  data  for  sample  Au-C6- 


Figure  24:  Low-frequency  noise  data  for  sample  Au-C6. 

than  500  Hz,  the  noise  is  white  noise.  At  frequencies  less  than  100  Hz,  the  noise  is  observed 
to  increase  with  increasing  bias  voltage. 

Figure  24  shows  the  corresponding  low  frequency  noise  (LFN)  as  a  function  of  fre¬ 
quency  for  the  various  bias  voltages.  The  LFN  was  obtained  by  appropriately  subtracting 
the  white  noise  from  the  measured  noise.  Figure  24  includes  the  power  fit  to  the  data 
using  Eq.  (5.4).  For  frequencies  greater  than  500  Hz,  the  noise  is  white  noise.  The  LFN 
can  be  fitted  by  1  / f1  with  the  exponent  7  =  0.80,  suggesting  that  with  decrease  of  the 
frequency,  the  LFN  increases  slower  than  a  typical  1//  noise,  i.e.,  the  LFN  is  suppressed 
as  compared  to  the  1//  noise. 
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6  Conclusions 


Our  theoretical  and  experimental  results  described  above  indicate  that  we  can  synthesize 
CLMPs  with  desired  particle  size  based  on  modeling  and  can  grow  high-quality  CLMP 
films.  The  TCR  of  the  CLMP  films  is  high  and  comparable  to  the  state-of-the-art  VO^. 
The  low-frequency  noise  of  CLMP  films  exchibits  a  weaker  frequency  dependence  than  the 
1  //  noise,  consistent  with  our  modeling  prediction.  This  weaker  frequency  dependence 
promises  a  low  noise  level  of  CLMP  films,  which,  together  with  the  large  TCR,  sug¬ 
gests  that  CLMP  films  have  a  great  potential  for  uncooled  microbolometers,  the  critical 
technology  for  our  soldiers. 
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